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An explicit perturbative solution to all orders is given for 
a general class of nonlinear differential equations. This so- 
lution is written as a sum indexed by rooted trees and uses 
the Green function of a linearization of the equations. The 
modifications due to the presence of zero-modes is considered. 
Possible divergence of the integrals can be avoided by using 
approximate Green functions. 



The Born expansion is a common tool of quantum 
mechanics. It states that, for a Hamiltonian H = 
Hq + V, a solution of Hip = Exp is given by (p — ipo + 
^^L (GF)™(/3o, where ipo is a solution of Hgipo = Etpo, 
and G = (E — Ho) -1 is the Green function correspond- 
ing to the unperturbed problem. The Born expansion has 
rendered innumerable services for pratical as well as the- 
oretical problems of quantum mechanics. A drawback of 
this expansion is that it is restricted to linear problems. 

In the present paper, the Born expansion will be ex- 
tended to non-linear problems. More precisely, we con- 
sider an equation of the type Fo(ip) = Fi(ip), where -Fb 
and F\ are functionals of tp. Fq describes the unper- 
turbed system and F\ its perturbation. For instance, the 
propagation of ideal optical solitons in an optical fiber is 
governed by the non-linear Schrodingcr equation, ppl 



.dtp ld 2 tp 2 



(1) 



We assume that we know a solution of the equation 
Fo(<Po) = 0. For Eq. ([!]), such a solution is the optical 
soliton 



?yexp (i(ry 2 — £ 2 )z/2 + i£r) 
cosh (t)(t - CO) 



where 77 and £ are parameters. 

In a true optical material, the ideal soliton is perturbed 
by various effects represented by the perturbation Q 

Fi(<p) = -il> + i(3-^ - at — + 

We will show that, if Fq and F\ are functionally dif- 
ferentiable, an explicit expression for the solution of 
Fo(ip(x)) = F\((p(x)) is the Butcher series 



ip(x) = <p (x) + V -4ttv(*; x). 

t a{t ) 



(2) 



The first part of the paper will be devoted to the ex- 
planation and derivation of Eq. (|^) , then the role of zero 
modes will be investigated. Finally, we show how ap- 
proximate Green functions can be used to avoid divergent 
integrals. 

The basic tool for the derivation of Eq.(|^) is that we 
sum over all rooted trees t. A tree is a graph without 
loop, and a rooted tree is a tree where one vertex is des- 
ignated as its root. We draw the root at the bottom of 
the tree. The smallest rooted tree is the single root •. 
The rooted trees with up to 4 vertices are 



I Y 



V Y 



Rooted trees have been introduced in 1857 by A. Cay- 
ley to represent derivatives of a function with respect to 
a parameter [||. In 1963, J.C. Butcher used rooted trees 
to write the solution of flow equations and derive new 
numerical methods Q. Since that work, series indexed 
by rooted trees are called B-series or Butcher series. In 
1998 Kreimer discovered that rooted trees underlie the 
basic mathematical structure of renormalization theory 

In Eq.(|2j), \t\ designates the number of edges oft. For 
instance, 1 1 1 = 2. To define the symmetry factor a(t), we 
need the collecting operator B + , ||] which starts from k 
trees t\, . . . ,tk and builds a new tree t = B+(t\, . . . ,tk) 
by joining the root of each of the k trees to a new ver- 
tex that becomes the root of t. For instance -B+(-) = I , 
B + ( • , • ) = V- The edges can be rotated around the ver- 
tices, so that 



.(.,!) = B + (l,.)=V 



Every rooted tree t (except the root) can be written as 
t = B + (ti, . . . , tk) for some ti,...,tk- Finally, a{t) is the 
symmetry factor of tree t. It is defined recursively by 

<K0 = 1. 



a(B + (C . . . ^)) = ni la(h) ni ■ • .n k \a(t k y 



(3) 



t^ fc ) means that t is obtained 



The notation t = B+tt™ 1 
by collecting n\ times tree ti,. . . , n k times tree where 
the k trees t\, . . . , tk are all different. 

The last term to define in Eq.(||) is (p(t;x). We recall 
that the functional derivative of F(ip(x)) with respect to 
a function ip(x) is M 

F(ip(x) + eijj{x)) - F(<p(x)) 



5F(<p(x)) 
5ip(x) 



lim 

e^0 



1 



This is called the Gateaux derivative in the mathematical 
literature JO. Then the functional derivative 

5F(<p(x)) . 6F(<p(x)) 

is defined as — — — — for ip(x) = o(y — x). 



My) 



8ip{x) 



From Fq and the unperturbed solution ipo (x) wc define 
the operator M with kernel M(x,y) 



M(x,y) 



SF (ip a (x)) 
6<p(y) 



Now we define the corresponding Green function G(x, y) 
by J dyM(x,y)G(y, z) = 6(x — z). The definition of 
G(x,y) involves the boundary conditions imposed on 
<p[x). A method to construct G(x,y) for solitons was 
proposed by Kawata and Sakai ^] . 

We can write Fo(<p(x)) as a sum over its functional 
derivatives: 



F (ip(x)) = F (ipo(x)) + J dyM(x,y){<p(y) - <po(y)) 
where 

6 n F (<p (x)) 



If 5 n l 

F 2 M,)) = ^-yd yi ...d y „^- 



) . ..6<p(y n ) 

x (<p(yi) - <Po(vi)) ■ ■ ■ (<p(Vn) ~ fo(Vn))- 

Since F (ip>o(x)) = 0, the equation F<o(<p) = Fi(<p) can be 
rewritten 



dyM{x,y)(ip(y) - <p (y)) - F{<p(x)), 



where F(ip(x)) = Fi(<p(x)) — F2(<p(x)). We operate the 
last equation by the Green function G to obtain 

<p(x) = (f (x) + J dyG(x,y)F((p(y)). 

This equation has exactly the form of Eq.(25) in Ref. 
p0| . and the proof given in Ref. |l(]] can be followed to 
show that <p(x) is given by Eq.(^|), where (p{t\ x) is defined 
recursively by 



ip(-;x)= J dyG{x,y)F((p Q (y)) 

for the root and, for t = B + (ti, . . . , ifc), by 

5 k F{Mv)) 



(p(t; x) = j dydzi . . . dz k G(x, y) 
xif(ti;z 1 )...(p(t k ;z k ). 



5ip(z!) . . . 5ip{z k ) 



To summarize, once the initial function ipo(x) and the 
corresponding Green function G(x, y) are known, the cal- 
culation of <p(x) up to any order is a mechanical appli- 
cation of a simple formula, which is well suited to com- 
puter algebra programs. To estimate the size of the terms 



in Eq.(||), we can use the fact that, if F\{(j)) is written 
as eF\ (</)), then (j>(t\x) is a sum of terms ranging from 

e (|t| + l)/2 t0 e \t\. 

To be complete, we must investigate the influence of 
zero modes. The term zero modes come from the theory 
of instantons in quantum field theory pTj , ^ . Within our 
framework, they are the solutions ip n (x) of the equation 
/ dyM(x,y)ip n (y) = 0. In other words, the zero modes 
are a basis of the kernel K of the operator M (x, z). 

In a linear problem, the zero modes are the solution of 
the unperturbed equation Hipo = Eipo , and the superpo- 
sition principle tells us that, if <p (rcsp. ip') is the solution 
given by the Born expansion starting from the zero mode 
ipo (rcsp. <p' ), then wc can start from ipo + (p' to obtain 
another solution (ip + ip'). The proper <po is determined 
by the boundary conditions. 

For a nonlinear problem, the superposition principles 
does not hold. For notational convenience, we assume 
that dim A" = 1 and -00 G K. We write the solution 
of the nonlinear problem as tp — (fo + p±_ + Xipoi where 
ip± G . As in the first part of the paper and in [fl0|| , 
we transform the original problem into an equation of the 
form 



<p(x) = <p (x) + F(ip(x)). 



(4) 



Then, such equation can be solved immediately using 
Butcher series, as was also noticed indepently by Schatz- 
man |l3| and Connes and Kreimer |Q. Notice that 
Butcher series can be obtained for quite general F(ip(x)): 
F can be a function of x (i.e. F{ip(x), x)), it can involve 
differentials, <p(x) and x can be multidimensional [ ^0[ . 

Now, we start from the equation Fo(ip) = Fi((p), and 
we propose a possible method to bring the problem into 
the form (Q). From ip = po + (p j_ + A^o we can write 

F (tp) = F a (tp Q ) + Mip±_ + XMtpo + F 2 {lp), 



Fi(yj) =F l (y ) + \- 



Sip 



where F 2 and F% are defined by the equations and 
5F 1 (p () (x))/Sip {x) = J dyip {y)SF 1 ((p (x))/5p(y) and 
M(x,y) =5F (<p Q (x))/5tp(y). 

Since Fq(pq) = and Mipo = 0, we obtain 



M<pj 



A 



SFxjcpo) 



F^tpo) + F(<p), 



where F = F3—F2. We need an independent equation for 
A, so we define the scalar product (/, g) = J dxf(x)*g(x) 
and we normalize ipo so that (ipo,ipo) = I. Then we use 
the fact that, for any g, (ipo,Mg) = to write 

-M^o,^|^) = Wo,F 1 (<p )) + (ih,F(<p)). 
Let a = —(ipo,5Fi((po)/5ipo), wc find 



2 



A = (i> ,F 1 (<p ))/a + {^ ,F{ v ))/a. 

To obtain the second equation, we define A as the pro- 
jector onto K- 1 , and we assume that A6F\((po)/Stpo = 0. 
Thus we obtain 

M(fj_ = A(Fi((p ) + F(<p)). 

If Go was the original Green function for M, we define 
G = AGoA and we obtain 

<Px = G(F 1 (( Po ) + F{<p)). 

We can group these two equations into a single one 

A = (tp ,Fi{<Po))/a + (4'o,F(ip + ip± + \ipo))/a 
tp ± = G{F l (<p )) + G(F(cp Q + <p ± + \iPo)). (5) 

This equation is now in the form required for the appli- 
cation of Butcher's method which writes the perturbative 
solution as 

t a w 

The zero-order terms are 

*S = (Vo,*i(¥>o))/a, 
& (x)= ( dyG(x,y)F 1 ( ( p (y))- 



For the roots, the functions < I> 1 and <£> 2 are defined by 
$ 2 (.;x)- / 6yG{x,y)F(My))- 



For a tree t = B+{t\, . . . , they are defined recursively 
by 



$\t)= J2 (V«) J ^Mv)^ 

01 — jk 

X& 1 (t 1 )...&>(t k ) 



S k F(<Po(v)) 



01— jk 

Xfc^ti)...*^). 

In the last expression, for ji = 1, then S/S& Ji = 8/Sipo(y) 

= $ 1 (*i)> and for 3i = 2 > thcn S/5&* = 
6/5<& 2 (zi), = 3> 2 (ij; Zi), and an integral is implic- 

itly assumed over the variable Zi- 

This was just an example of the general strategy avail- 
able for the treatment of zero modes. In specific prob- 
lems, it might be more efficient to take also account of 
the dependence of F ((p) on A. 



As a last point, we want to show on a simple exam- 
ple, that using a modified Green function can make the 
problem much better behaved. The example we want to 
discuss is 



d 2 (p(x) 
dx 2 



+ <p(x) - 2gip(x). 



This problem as a nontrivial solution 

, x) = J L_ 

yfg cosh x ' 

The linearization gives 

= ( _^ + 1 _ 



8tp(y) 



— 72-)S(x-y). 
cosn x 



The equation Mtp = has two solutions: a normalizablc 
solution ipo(x) and a non normalizable one ^i(x). 



tp (x) 



sinha; 
cosh 2 x 



. . x sinh x — cosh x 

ipi(x) = 3 g 1" cosh a;. 

cosh x 



The kernel of M is the one dimensional subspace gen- 
erated by tpo(x). By Wronski's method ( @, p. 900), 
the Green function for M is 

G°{x,y) = -^V'o(a;<)'0i(^>), 

where, (x < ,x > ) = (x, y) if x < y, and (a5 <( ar>) = (y,x) 
if a; > y. However, this Green function is not well 
behaved because of the term cosh a; which gives diver- 
gent integrals. If calculations are made with G°(x,y), 
the apparatus of renormalization theory has to be used. 
This can be avoided by defining a well-behaved function 
G 1 = G° + 5G, where 8G(x,y) = M%<) cosh(x>)/2. 
It can be checked that G 1 {x,y) is now exponentially de- 
creasing. The equation AG M(f± = ip± becomes 

AG 1 Mip ± =ip±- A6GMip±. 

Therefore, Eq.(|J) becomes 

<p ± = AG 1 A(F 1 (ipo) + F(ip Q + <P± + A^o)) 
-ASGMA(pj_- 

Note that, in the sense of distributions 

3 



dzSG(x,z)M(z,y) 



-1 



cosh 



2 cosh 
-6{x-y) 



5{x - y) 



is no longer divergent when integrated over x. 

A method was proposed to write the general term of 
the perturbative solution of a nonlinear problem. This 



3 



method is a nonlinear generalization of the Born expan- 
sion. Although Butcher scries arc not widely known in 
Physics, they are commonly used in Numerical Analysis, 
where they have proved their power fln^ , |i~7I] , and where 
many results are now available. 

From the theoretical point of view, Butcher series en- 
able us to write the general term of the perturbative so- 
lution of nonlinear equations. And since we can write it, 
we can manipulate it to investigate stability, long-time 
behavior, convergence, bifurcation, transition to chaos, 
resummation, and so on. 

From the practical point of view, the perturbation of a 
soliton is usually calculated by solving recursive differen- 
tial equations, and the computation becomes very cum- 
bersome beyond the first order 0]. By contrast, Butcher 
series offer a systematic method which is easily imple- 
mented on a computer. 

In Ref. [S, several perturbations of the optical soli- 
ton are reviewed. All can be considered as induced by 
a change of variable. Such a change of variable gives a 
modified nonlinear equation with a modified linearization 
M' . The Butcher series solution gives the general term 
of the perturbative solution for this modified problem. 
Moreover, the non-linear change of variable can be also 
implemented with trees. Jl8| This, together with the fact 
that rooted trees have recently appeared in singularity 
theory [|l9j, renormalization of quantum fields [^0|, and 
non-commutative geometry plj indicate that they pro- 
vide us with a tool particularly well adapted to nonlinear 
problems. 
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